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Abstract 

Meniscus  oscillation  usually  occurs  after  jetting  drops  in  Drop-on-Demand  inkjet  head.  It  is  important  for  ink  refill  motion. 
Ink  refill  motion  affects  many  jetting  performance  in  inkjet  technology.  Therefore  it  is  very  important  to  understand  ink 
refill  motion  process  for  designing  inkjet  head  dimensions  and  ink  properties,  such  as  nozzle  diameter,  barrier  thickness,  ink 
viscosity  and  surface  tension.  This  meniscus  oscillation  is  a kind  of  under-damped  oscillation.  The  refill  time  is  defined,  as 
the  time  required  returning  to  the  initialization,  the  free  surface  obeys  damped  oscillation  and  oscillates  between  meniscus 
mounding  and  recession  when  the  amplitude  of  the  oscillation  reduces  to  dispersion.  This  study  researches  several  inkjet 
heads  using  computational  fluid  dynamics  simulation.  CFD  solver,  such  as  FLOW  3D,  is  used  to  solve  the  problem.  After 
calculation,  the  results  are  plotted  with  post-processor  such  as  plot  software  and  output  to  printers.  This  paper  shows  the 
break-off  time,  the  refill  time  and  operation  frequency  of  the  inkjet  head. 

Keywords : Thermal  inkjet,  Meniscus  oscillation,  CFD,  fluid  dynamics. 

1*  Introduction 

In  recent  years,  an  inkjet  printer  has  become  the  most  dominant  output  device  of  PC  due  to  its  low  cost,  high  print  quality 
and  color  printing.  However,  the  print  quality  is  closely  related  to  the  ink  droplet  generation  process.  The  ink  droplets  are 
usually  jetted  by  a pressure  source  generated  by  thermal  bubbles  or  piezoelectric  transducers.  When  pressure  is  applied  to 
the  bottom  of  the  chamber,  the  liquid  inside  the  diamber  is  compressed  in  order  to  form  droplets  jetted  from  the  nozzle 
which  is  on  the  top  of  the  chamber.  After  the  droplet  breaks  off  at  the  nozzle,  the  chamber  is  then  refilled  with  ink  again. 

The  refill  motion  of  ink  affects  many  printing  performances.  In  1977,  Beasely1  calculated  the  refill  motion  for  a 
piezoelectric-driven  inkjet  head  based  upon  a simple  mathematical  model,  and  compared  with  two  experimental  results.  The 
analysis  reveals  that  the  equivalent  length  is  much  shorter  than  the  physical  length  during  the  refill,  and  the  calculated  refill 
times  agree  with  the  experiments  observed.  1981,  Kyser2  introduced  a new  mathematical  model  which  can  be  applied  to  the 
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known  embodiments  of  jets  and  shown  meniscus  oscillation  in  the  refill  motion.  Saitor3  used  a simplified  model  under  some 
very  simple  assumptions  to  study  the  refill  time  and  oscillation  characteristics.  In  his  study,  the  neglect  of  the  drop  ejection 
process  caused  some  disagreements  with  experimental  results. 

This  study  uses  a Computational  Fluid  Dynamics  (CFD)  model  to  simulate  the  drop  jetting  behave  of  an  inkjet  head,  and  to 
observe  the  meniscus  oscillation  of  free  surface  near  the  exit  of  the  nozzle  plate  after  drop  jetting. 


2.  Numerical  description 

Simulation  of  an  inkjet  head  is  one  of  the  challenging  tasks  in  CFD.  The  structure  size  is  quite  small,  in  the  order  of  10 
micrometer,  and  causes  severe  truncation  errors  in  the  calculated  results.  Another  challenge  is  keeping  track  of  the  free 
surface  between  liquid- gas  phases.  Free  surface  plays  a very  important  role  in  our  simulations  because  modeling  requires 
tracking  the  interface  in  a more  accurate  resolution.  CFD  to  be  used  in  this  case  becomes  very  complicated,  the  code  has  to 
be  capable  of  capturing  the  interface  accurately.  Some  packages  are  capable  of  dealing  the  density  discontinuity  of  the  free 
surface  by  smoothing  the  interfaces  over  a few  grid  cells,  but  not  accounting  for  the  sharp  changes.  In  order  to  simulate  the 
meniscus  oscillation  of  the  inkjet  head,  the  code  has  to  have  the  capabilities  of  tracking  free  surfaces,  keeping  its  sharp 
changes  in  one  cell  distance  and  evaluating  the  free  surface  curvatures. 

VOF4  (Volume  Of  Fluid)  method  utilizes  a finite  difference  to  represent  the  free  surfaces  and  interfaces  which  are  arbitrarily 
oriented  with  respect  to  the  computational  grids.  Flow_3D  is  a CFD  solver  using  the  VOF  method  in  two-phase  flow.  It  can 
track  free  surfaces  very  accurately,  and  has  a certain  degree  of  accuracy  in  surface  sharp.  The  VOF  method  is  a wide-use 
two-phase  flow  model  to  tackle  interface  between  gas  and  liquid  phases.  VOF  has  been  implanted  into  many  packages  such 
as  SOLA,  SOLA- VOF,  NAS  A- VOF  and  so  on.  Today,  many  inkjet  manufacturers  and  developers  use  those  packages  to 
simulate  the  jetting  process. 

VOF  has  defined  a function  F whose  value  of  unity  at  any  point  occupied  by  fluid  and  zero  otherwise.  The  average  value  of 
F in  a cell  then  represents  the  fractional  volume  of  the  cell  which  is  full  of  fluid,  while  a zero  value  indicates  that  the  cell 
contains  no  fluid.  There  exists  a free  surface  while  the  value  of  F is  between  zero  and  one  in  a cell.  Although  VOF  can 
locate  the  free  boundary  nearly  as  well  as  a distribution  of  marker  particle  method4,  and  has  the  advantage  of  a minimum  of 
stored  information,  the  method  is  worthless  unless  an  algorithm  can  be  devised  for  accurately  computing  the  evolution  of 
the  F field.  The  time  dependence  of  F is  governed  by, 


8F  8F  dF  „ 

h u h v — = 0 

dt  dx  dy 


(1) 


Eqn.l  describes  that  F moves  with  the  fluid,  and  is  the  partial  differential  equation  of  marker  particles.  Thus,  the  VOF 
method  provides  a simple  way  to  track  the  free  surface. 


3.  Results  and  Discussions 

The  siiriulation  of  Inkjet  head  has  some  difficulties  such  as  the  criterion  of  the  bubble  formation  under  the  condition  of  super 
heating  and  the  free  surface  capturing.  The  treatment  of  the  thermal  bubble  growth  of  an  inkjet  head  belongs  to  the 
super-heating  model.  So  far,  there  are  no  commercial  packages  can  deal  with  super-heating.  In  fact,  it  is  a big  challenge  for 
the  inkjet  simulation.  This  study  introduces  a solid  piston  to  substitute  for  the  thermal  bubble,  which  is  a pressure  source  for 
inkjet  jetting  process.  Fig.l  is  the  contour  of  computational  domain  for  Flow_3D.  The  computational  domain  is  an 
two-dimensional,  axisymmetric  inkjet  head  which  the  radius  of  the  nozzle  plate  is  30  fim  and  50//m  in  thickness.  The 
height  of  barrier  (or  ink  supply  channel)  is  30  fi  m.  For  the  convenience  of  calculation,  the  solid  bottom  including  a movable 
piston  has  been  set  as  2 fim  thickness.  The  void  region  is  500 //m  in  length,  52 /m  in  radius.  The  piston  has  a specified 
velocity  as  a function  of  time  shown  in  Fig.2.  The  boundary  conditions  are  given  as  follows:  symmetry  in  central  line,  wall 
condition  in  bottom,  continuative  condition  on  the  top  exit  and  one  atmosphere  pressure  condition  on  the  right.  The 
viscosity  coefficient  and  surface  tension  coefficient  are  3.5  cps  and  50  dyne/cm2,  respectively 

Fig.  3 shows  the  variation  of  the  interface  height  in  the  axial  direction  for  the  case  in  which  the  nozzle  radius  is  30  fi  m and 
the  height  of  barrier  is  30  fi  m.  The  lowest  point  of  27  fi  sec  is  the  break-off  time  of  the  droplet.  This  figure  shows  that  the 
oscillation  of  the  free  surface.  Table  1 shows  the  computational  results  for  four  kinds  of  nozzle  radii  and  four  kinds  of 
barrier  heights.  From  the  results,  we  found  that  the  smaller  the  nozzle  radius,  the  less  the  break-off  time.  Besides,  the  barrier 
height  does  not  influence  the  break-off  time  if  the  nozzle  radius  is  kept  constant.  The  results  of  the  refill  time  show  that  the 
larger  the  nozzle  radius,  the  longer  the  time. 

Some  of  the  results  do  not  obey  this  tendency,  we  found  that  small  satellite  droplets  separated  from  main  droplet  or  second 
droplet  collide  back  to  the  free  boundary  after  jetting.  This  may  cause  slightly  large  oscillation  on  the  free  boundary  and 
destroy  the  natural  frequency  of  damped  oscillation.  We  observe  from  experiments  that  there  are  usually  a small  amount  of 
ink  deposits  around  the  muzzle.  Perhaps,  the  residual  ink  around  the  muzzle  may  have  the  similar  phenomena  compared 
with  the  simulation  results.  The  reason  is  that  the  smaller  droplet  has  opposite  momentum  while  separating  from  the  larger 
droplet  and  has  larger  buoyancy  effect. 

Table  1 also  shows  the  oscillation  frequency  of  the  results.  Almost  the  frequencies  are  larger  than  2kHz  and  smaller  than 
3kHz.  This  shows  the  structure  can  work  at  3kHz. 

From  the  above  results,  we  use  several  kinds  of  nozzle  radii  and  barrier  heights  to  calculate  the  break-off  time  of  the  droplet 
and  the  ink  refill  time.  The  results  showr  that  CFD  can  use  to  simulate  the  behaviors  of  inkjet  heads.  Although  there  are  some 
differences  in  the  break-off  time  and  the  refill  time,  the  operation  frequencies  of  the  inkjet  heads  still  work  about  2kHz. 
Here,  the  design  operation  frequency  is  3kHz. 
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Fig.l  2D  computational  domain 
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